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What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum 
system? We provide an efficient algorithm to simulate any desired two-body Hamiltonian evolution 
using any fixed two-body entangling n-qubit Hamiltonian and local unitaries. It follows that uni- 
versal quantum computation can be performed using any entangling interaction and local unitary 
operations. 
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A central goal of quantum physics is to understand and 
control quantum dynamics. Recently, the emergence of 
fields such as quantum control , laser cooling 
turn communication and quantum computation 
focused efforts to understand and control quantum dy- 
namics at the single quantum level. 

Our interest is in the dynamics of composite quantum 
systems. An especially important example of such a sys- 
tem is a quantum computer, which is a composite of a 
large number of two-level quantum systems (qubits). We 
wish to determine which interactions are sufficient for the 
simulation of arbitrary quantum dynamics in such a sys- 
tem. Our results demonstrate equivalence between this 
property of universality and the ability to entangle all 
components of the system. 

More precisely, we consider the following problem: 
what dynamics can we produce with a specified two- 
body, n-qubit Hamiltonian, given the ability to perform 
arbitrary local unitary operations on individual qubits? 
Under these conditions, we exhibit an explicit algorithm 
which shows that any Hamiltonian which produces entan- 
glement can be used to efficiently simulate an arbitrary 
two-body dynamical operation. This holds even if the 
Hamiltonian alone is only capable of producing a small 
amount of entanglement. 

It follows that any entangling interaction, together 
with local unitaries, is sufficient to perform universal 
quantum computation. Our results thus confirm the folk- 
lore belief that it is the ability to entangle that is the 
crucial element in quantum computation. 

Substantial prior work has been done on universal op- 
erations, and many specific sets of universal gates are 
known |,|. Our work differs from previous work on the 
general requirements for universality in several regards. 
Closest is the work in || and j?]] , where it was shown that 
almost any two-qubit quantum gate is universal for quan- 
tum computation. This work focused on unitary gates 
rather than continuous-time Hamiltonian evolution, and 
did not explicitly determine which sets of unitary gates 
are universal. Our work explicitly determines which two- 



body Hamiltonians, together with the additional require- 
ment of local unitaries, are universal. Furthermore, in Q 
and it was assumed that gates could be independently 
applied to any pair of qubits in the computer, and thus 
required the ability to turn on and turn off interactions 
between different pairs of qubits. By contrast, we assume 
only a fixed entangling operation. 

Our techniques make use of generalisations of stan- 
dard nuclear magnetic resonance (NMR) techniques for 
decoupling and refocusing ||||. Similar ideas have been 
applied by (H]|n|] to the problem of efficiently implement- 
ing coupled logic gates using a restricted class of hamil- 
tonians which arises naturally in NMR. 

The structure of this Letter is as follows. We begin 
with a precise formulation of our goals and results. A 
specific two-qubit example is given to illustrate our tech- 
niques, and the general algorithm is described for the 
case of an arbitrary two-qubit system. The efficiency of 
the algorithm and the effect of errors are then discussed. 
We conclude by generalizing the algorithm to n-qubit 
systems. 

An arbitrary Hamiltonian on n qubits can be given the 
operator expansion 



H - ^2 h n---u a n 



'3n> 



(1) 



where the hi 



are real numbers and o\ , <y% , 03 are the 



usual Pauli sigma matrices, with ctq = I the identity. Our 
discussion is restricted to the case of time-independent 
Hamiltonians containing only one- and two-body terms, 
that is, if hj 1 _,j n ^ then only one or two of the ji, . . . , j n 
are not equal to zero. If the Hamiltonian contains a non- 
zero contribution to <8 07 then we say the Hamiltonian 
couples systems k and I. This focus on two-body Hamil- 
tonians is a mild restriction as most candidate systems 
for quantum information processing are of this type. 

Under what circumstances is it possible to produce en- 
tanglement between an arbitrary pair of systems, even 
ones that are not directly coupled by the Hamiltonian 
HI Not surprisingly, Hamiltonians which have terms 
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coupling systems k and I can produce entanglement be- 
tween these systems. We say that systems k and k' are 
connected if there is a sequence (k, k\, . ■ ■ , k m , k') such 
that each adjacent pair in the sequence is coupled by H. 
It is clear that if k and k' are not connected then no en- 
tanglement can be created between them, and thus it is 
not possible to perform an arbitrary unitary operation on 
the system. Conversely, it follows from our later discus- 
sion (and is a priori plausible) that if a pair of systems 
is connected then it is possible to create entanglement 
between them (cf. [jl2|,|l3|] ) . This motivates our defini- 
tion of a two-body entangling Hamiltonian as a two-body 
Hamiltonian such that all pairs of systems are connected. 
The main result of this Letter is the following: 

Let H be a given two-body entangling Hamilto- 
nian on n qubits, and let if be a desired two-body 
Hamiltonian on n qubits. Then we have an efficient 
algorithm to simulate evolution due to K using only 
(a) the ability to evolve according to H and (b) the 
ability to perform local unitaries on the individual 
qubits. 

In particular, given such a Hamiltonian it follows that we 
can perform an arbitrary two-qubit unitary gate on any 
specified pair of qubits. Thus, by well-known universal- 
ity results , we may efficiently perform any quantum 
computation. 

Three elementary observations about Hamiltonian evo- 
lution form the key to our methods: 

(A) Imagine we can evolve according to the Hamiltonian 
J, and perform unitary operations U and UK Then it fol- 
lows from the identity e - ltuju ^ — JJe~ %tJ W that we can 
exactly simulate evolution according to the Hamiltonian 
UJW. 

(B) Imagine we can evolve according to Hamiltonians J\ 
and J2. Then we can simulate evolution due to J\ + J 2 
for small times A, due to the approximate identity 

e -iA(Ji + J 2 ) e -iAJ l£ -iAJ 2 /g\ 

Initially we treat this identity as though it is exact, and 
analyze the effect of errors later. 

(C) Imagine we can evolve according to a Hamiltonian 
J. Then, by appropriate timing, we can exactly simulate 
evolution according to AJ for any A > 0. 

The basic idea can be illustrated using a two-qubit ex- 
ample. Suppose we have the ability to evolve according 
to the two-qubit Hamiltonian: 

H = Z (gi I + 2X ® Z + Z <g> Z, (3) 

where X, Y and Z are a convenient shorthand for the 
Pauli sigma matrices. H couples the two qubits, and is 
thus a two-body entangling Hamiltonian. The first step 
of our procedure is to show that H and local unitaries 
can be used to simulate evolution according to the largest 



coupling term in H, in this case X ® Z. This follows im- 
mediately from the identity 

X®Z= ~(X(g>I) J H'(X(g>I)t + 2#, (4) 

and our earlier observations (A), (B) and (C). 

Using the ability to simulate evolution by the Hamil- 
tonian X (& Z we can easily obtain the ability to simulate 
a Hamiltonian which is any product of Pauli matrices. 
Products of the form I ® (7j and o~j ® I follow immedi- 
ately from our ability to do local unitaries. Products of 
the form Oj <S> cr k follow from observation (A) and the fact 
that o j ®o k = {U® V)X ® Z(U ® V)* for appropriate 
single-qubit rotations U and V. It is easy to see that 
observation (A) also allows us to simulate terms of the 
form — o~j <® <7fc. An arbitrary two-qubit Hamiltonian K 
can be decomposed as a linear combination of products 
of Pauli matrices, and thus by observations (B) and (C) 
may be simulated using our ability to simulate X ® Z . 
Thus K may be simulated using H and local unitaries. 

The general two-qubit case follows using similar tech- 
niques. Suppose H is an entangling Hamiltonian and 
choose r,s / such that \h r>s \ is maximized. It is easy 
to verify that 

(o-j ® o- k )H(o-j ® a k y 

j€{0,r},*e{0,*} 1 ' 1 

hrfiVr ® I + ho, s I ® OS + h . Q I <g> I 

Using observations (A), (B) and (C) it follows that both 
ay ® cr s and — a r ® a s can be simulated using H and 
local unitaries, and thus the result follows for a general 
two-qubit Hamiltonian K. 

In more detail, suppose we wish to simulate K for a 
non-infinitesimal time t > 0. We have shown that we 
can approximate evolution due to K for a small time A 
by applying an appropriate sequence of evolutions due 
to H and local unitaries. Such a simulation requires, 
in general, 36 separate periods of evolution due to H, 
interleaved by single-qubit unitary gates applied to the 
two qubits. To simulate the evolution due to K over a 
time t we break the interval t into N increments of length 
A = t/N, and perform the simulation of K for each incre- 
ment, repeating the small time-step procedure N times, 
for a total of at most 36iV separate periods of evolution 
due to H. 

Let us turn to the sources of error inherent in our sim- 
ulation procedure. The procedure uses only the observa- 
tions (A), (B) and (C). (A) and (C) are in principle exact, 
however the identity Eq. (||) used in (B) only holds ap- 
proximately. In order to perform a good simulation of 
K we therefore need to choose a timestep A sufficiently 
small that Eq. (Q) is a good approximation. 

To do the error analysis, we introduce a measure quan- 
tifying how well our simulated evolution approximates 
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the desired evolution due to K. That is, we wish to com- 
pare the unitary evolution W' achieved by our simulation 
with the unitary evolution W — exp(-iKt) that we wish 
to simulate. We use as our measure of error the operator 
norm of the difference between W and W, \\W — W'\\, 
defined by \\A\\ = max^.uui This is physically 

well-motivated since two operators W and W are close 
according to this norm if and only if the difference in 
their effects on an arbitrary state is bounded by a small 
number. The actual measure of error used is not all that 
important, but we find it useful to demand the following 
two properties, both of which are satisfied by the opera- 
tor norm: (1) stability under tensor product with ancilla 
systems, that is, \\A\\ = \\I®A\\; and (2) invariance under 
unitary transformations, that is, ||A|| = ||V'AV|| for any 
unitaries V, V . This latter property implies the chaining 
inequality for any unitary operators V\, V2, W\, W%: 

|| Wi - V 2 W 2 \\ < || Vi - V2W + \\Wi - W 2 ||- (6) 

We bound the errors induced by the approximation in 
Eq. (||) using the inequality Q : 

|| e -ir(Ai+...+A m ) _ e ~ iTA i _ _ _ e ~irA m || 
l<j<fe<m 

where r is a positive real number and the Aj are Her- 
mitian operators. Applying this bound and the chaining 
property to the procedure we've described gives 

\\W -W\\ <CD 2 tA, (8) 

where C is a constant which we can easily bound to 
be at most 10 4 , and D is a parameter determined 
by the properties of H and K as follows. Let h = 
maxij \hij\,k = m&Xij \kij\, where K has the operator 
expansion . fejjcrj ® oj. Then D = \hk/h rtS \. 

The error bound Eq. (||) can be improved substantially 
in several ways. The linear dependence on A in Eq. (|^) 
is due to the technique used to simulate sums of Hamil- 
tonians, namely e~ lA ( Jl+j2 ) = e -*A Ji e -iAJ 2 + q(A 2 ). 
Each simulation step thus contributes an error 0(A 2 ), 
and there are t/A such steps for a total error 0(tA). 
Higher-order approximation techniques ]i~4| ] can be used 
to obtain more accurate simulations. For example, iden- 
tities such as 

e -iA(,/ 1 + J 2 ) _ e -iA.h/2 e -iA.h e -iAJ 1 /2 + Qr/^3\ ^ 

yield a cumulative error which is 0(iA 2 ). In general, an 
approximation analogous to Eq. (^) but accurate to order 
A fc leads to a cumulative error 0(iA fe_1 ). The tradeoff 
is such that higher-order approximations require the use 
of somewhat more complicated gate sequences for each 
timestep. In practical applications, this additional com- 
plication must be balanced against the improvement in 
accuracy to achieve optimal results. 



A second way to improve the bound in Eq. (||) is 
to leverage specific knowledge of the given and desired 
Hamiltonians. For example, imagine that we have avail- 
able the Hamiltonian of Eq. (||), and wish to simulate a 
controlled-NOT gate ||. We can do this more efficiently 
than implied by the identity in Eq. (|5|) by examining the 
properties of the controlled-NOT. Up to an unimportant 
global phase, the controlled-NOT may be generated by 
applying the Hamiltonian / <g>A + Z®/ — Z ® A for a 
time t — 7r/4. The terms in this Hamiltonian commute, 
so we have controlled-NOT = e -i(i^x)t e -i(zm)t e i(z«,x)t_ 

Thus, to simulate the controlled-NOT for a time t, it suf- 
fices to simulate evolution according to the Hamiltonian 
K = —Z ® X, followed by local unitaries. We observe 
that 

K=(R®RX)( X ® I W X 4 ®V +H (R®RX)i (10) 

where R is the Hadamard gate ||, denoted here by R 
instead of the usual H to avoid confusion with the given 
Hamiltonian. Using the method outlined earlier gives a 
cumulative error 8tA. If we wish to have an accuracy 
of 10 -3 this corresponds to roughly 10 4 periods of evo- 
lution according to H, interleaved with local unitaries. 
This number of operations is probably too large to be 
practical, however it is substantially better than is ob- 
tained using the general bound Eq. (||). 

Further improvement may be obtained by using the 
higher-order approximation Eq. (|^). Using the op- 
erator norm, simple algebra shows that the correc- 
tion in Eq. (^|) may be bounded to order A 3 by 
lll^illll J2W (|| Ji|| +2||J 2 ||)A 3 . In this specific example, 
this reduces to ||ff|| 3 A 3 for a cumulative error of at 
most -^\\H\\HA 2 . Bounding \\H\\ by \\Z <gi 7|| + 2\\X O 
Z|| + ||Z <g> Z\ =4 we see that the cumulative error is 
at most \tA 2 . Therefore, to achieve an accuracy of 10 -3 
in our simulation of the controlled-NOT we need approx- 
imately 10 2 periods of evolution due to H, interleaved 
with local unitaries. Further improvements may be ob- 
tained by using better approximations than Eq. ([)]) . 

The number of operations required to simulate an ar- 
bitrary unitary operation can thus be substantial. In 
practice, this disadvantage may be offset by the advan- 
tages gained in using the natural coherent interactions 
present in a system. Furthermore, our results merely 
provide a lower bound on the efficiency with which it is 
possible to simulate an arbitrary unitary operation, and 
provide substantial impetus to search for better methods 
in specific cases. 

We now turn to the n-qubit case. The basic idea is to 
reduce the problem to the two-qubit case already solved. 
We divide the system into two parts, a principal system P 
consisting of two qubits which are coupled by the Hamil- 
tonian H, and the remainder of the system, denoted S. 
We use a technique generalizing the work in |0,fq| that 



3 



turns off all interactions between P and S and within S, 
leaving only the interactions present in P. These inter- 
actions can then be used, as before, to simulate arbitrary 
dynamics on the two qubits in P. Thus it is possible to 
simulate arbitrary dynamics on any two qubits coupled 
by the Hamiltonian H . Finally, an arbitrary interaction 
between qubits k and k' may be effected by performing 
a sequence of SWAP gates between the qubits connect- 
ing k and k' , applying the desired interaction, and then 
swapping back. 

The first step is to decouple systems P and S. To do 
this, let Xs denote a tensor product of X operators ap- 
plied bitwise to all the qubits in S. Define Y$ and Z$ 
similarly. Observe that forming the Hamiltonian 



H + XsHXl + YsHY^ + ZsHZl] (11 



leaves the Hamiltonian on P invariant, but eliminates all 
coupling terms between P and S, and all single-system 
terms on S. 

We now explain a recursive construction to eliminate 
all remaining couplings in S. First, we break the block 
S into two blocks So and Si of approximately equal size. 
We decouple So and Si by forming the Hamiltonian 



H": 



H' + X Sa H'X 



L S 



Y So H'Yl + Z So H'zl 



(12) 



Next, we break So into two blocks Soo and Soi of approx- 
imately equal size, and break Si into two blocks Sio and 
Sn of approximately equal size. We can decouple Soo 
from Sio, and Soi from Sn in a single step by forming 
the Hamiltonian 



H'" = - 

4 



H" 



x S 00 X SlQ H" x\ aa X 



+ Y Sno Y Slo H"YlYl > + Z Sno Z Slo H"zl oQ Zl ic ] . (13) 

We repeat this blocking procedure \log 2 (n — 2)] times 
to decouple all the terms in S, leaving a sum over 
O(4 log2 ") = 0(n 2 ) terms involving the conjugation of 
H by local unitary operations. 

Thus, simulating a Hamiltonian K applied to P for a 
time t requires the use of 0(n 2 ) periods of evolution due 
to H , interleaved with local unitaries. Using a similar 
error analysis to that described earlier, and the stability 
property of the operator norm, we find an error 0{n 2 tA). 
In practice it may be possible to do substantially better 
by leveraging our knowledge of specific systems, and us- 
ing better approximations. 

A number of problems will be addressed in future work, 
including: (a) the extension of our results beyond the 
qubit model to higher-dimensional systems; (b) the fault- 
tolerance Jl6|,p| of our simulation techniques; (c) the op- 
timization of our techniques for specific systems; (d) the 
further study of the general requirements for universal 



computation (cf. Jl7|~p0t). For example, it is likely inter- 
esting to impose restrictions on the class of local unitary 
operations that may be applied during the computation, 
perhaps adopting a cellular automata model in which op- 
erations are applied nearly homogeneously across the en- 
tire system. 

The results presented in this paper demonstrate that 
all two-body, n-qubit entangling Hamiltonians are equiv- 
alent in the sense that any such Hamiltonian can be used 
to efficiently simulate any other with the aid of local uni- 
tary operations. We conjecture that the same result is 
not true for /c-body Hamiltonians where k > 2. In the 
case k = 3, for example, it is possible to construct a 
Hamiltonian which, alone, can only generate GHZ-type 
entanglement. It would be surprising if this Hamiltonian 
could be used to simulate a Hamiltonian which produces 
W-type entanglement because these types of entangle- 
ment appear to be fundamentally different ]2l| ] . It would 
be of interest to determine, in general, what characteris- 
tics of two sets of Hamiltonians determine whether they 
are equivalent. 

We thank Carl Caves, Andrew Childs, Ike Chuang, 
Chris Dawson, Ivan Deutsch, and Tobias Osborne for 
helpful discussions. 
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